Nonlinear mechanics with photonic crystal nanomembranes 
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Optomechanical systems close to their quantum ground state (QGS) ij, @] and nonlinear nano- 
electromechanical systems (NEMS) are two hot topics of current physics research. Demonstrating 
the QGS allows to shed new light on quantum coherent effets in meso- and macroscopic systems, 
whereas NEMS operated in a nonlinear regime [3] are used either to demonstrate the underlying 
physics of sheer nonlinear effects [J-[ylor as RF amplifiers, with sensitivity improvement both at the 
classical 0,0] or quantum noise level [3, [13]. As high- reflectivity and low mass are crucial features to 
improve optomechanical coupling towards the QGS, we have designed, fabricated and characterized 
photonic crystal (PhC) nanomembranes at the crossroad of both topics. Here we demonstrate 
a number of nonlinear effects with these membranes. We first characterize the nonlinear behavior of 
a single mechanical mode and we demonstrate its nonlocal character by monitoring the subsequent 
actuation-related frequency shift of a different mode. We then proceed to study the underlying 
nonlinear dynamics, both by monitoring the phase-space trajectory of the free resonator and by 
characterizing the mechanical response in presence of a strong pump excitation. We observe in par- 
ticular the frequency evolution during a ring-down oscillation decay, and the emergence of a phase 
conjugate mechanical response to a weaker probe actuation. Our results are crucial to understand 
the full nonlinear features of the PhC membranes, and possibly to look for nonlinear signatures of 
the quantum dynamics [12i] . 



We focus on the mechanical nonlinear behavior of pho- 
tonic crystal nanomembranes developed as end mirrors 
in a Fabry-Perot optomechanical cavity [11]. Details 
about the mechanical and optical properties, as well 
as the fabrication process of these membranes are de- 
scribed in previous publications [ll|, [isl, [l^. They are 
10 /im X 20 /im x 260 nm indium phosphide slabs sus- 
pended over the substrate by four decoupling bridges 
(0.5-/im width and 6- to 12-/im lengths) located at the 
nodes of the mechanical mode of interest [l5|, 5.9 /im 
away from the center of the membrane. 

To probe the mechanical response of the membrane, 
the sample is actuated by a piezoelectric stack, and the 
membrane displacement is monitored by a Michelson in- 
terferometer (see Fig. [T|). The sample is operated in a 
vacuum chamber at 10~^ mbar to prevent from air vis- 
cous damping and squeezed film effects. We use three 
different membranes with resonance frequencies Qo/27r = 
782 kHz, 1036 kHz, and 1057 kHz, and quality factors 
Q 000, mainly limited by surface effects. 

A network analyzer actuates the sample and moni- 
tors the resulting displacement spectrum. Figure [2] (left) 
shows the forced oscillation amplitude Xp[l]p] as a func- 
tion of the actuation frequency f^p, for different actua- 
tion powers. The mode of interest, close to 1 MHz for the 
membrane used here, exhibits a strong nonlinearity even 
for an actuation voltage limited to a few volts applied to 
the piezoelectric stack, corresponding to displacements 
of a few picometers only for the substrate and at the 
nanometric level for the membrane. 
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FIG. 1. Interferometric setup: we use either a network ana- 
lyzer to probe the mechanical response of the membrane (a)^ 
or a signal generator and a spectrum analyzer triggered by 
the generator to test the ring-down behavior (h)^ or both of 
them to perform pump-probe experiments. The inset shows 
a typical photonic crystal membrane. 



The nonlinearity can be accounted for by the Dufiing 
model in which the resonance frequency has a quadratic 
dependence with the displacement x^{t): 

Xp{t) + Tx^it) + [1 + f3xl{t)] x^{t) = a^{t), (1) 

where a^{t) = 2apCos(l]pt) is the driving force, F = 
l^o/Q the mechanical damping, and P the nonlinearity 
strength. 

At first order the forced displacement is mainly 
monochromatic at the actuation frequency l^p, with a 
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FIG. 2. Oscillation amplitude Xp[Qp] of the membrane as a 
function of the piezoelectric actuation frequency Qp. Linear 
to nonlinear transition is clearly visible on left curves a to e, 
obtained with an upward frequency sweep of the network ana- 
lyzer and for increasing actuation powers (10 dBm to 30 dBm 
with a 5-dB step). Curves on the right are swept either up 
or down with the frequency generator at 30 dBm and exhibit 
a typical hysteresis cycle as expected from the theoretical fit 
(dashed curve). 



Fourier component given by (assuming Q ^ 1) 



nl{i + e)-ni-iTnQ' 



(2) 



This expression corresponds to the usual lorentzian reso- 
nance, except for a shift of the resonance frequency pro- 
portional to the mean square displacement through the 
term 



(3) 



This shift is responsible for the spectra observed in Fig. 
El as shown by the theoretical fit on the right (dashed 
curve). Although higher order effects must be consid- 
ered to completely fit the displacement, values of {3 can 
nonetheless be evaluated in the range of 10^^ m~^. Also 
note that such a more complete derivation [l6| would 
lead to the emergence of odd harmonics of motion with 
an amplitude at least e ^ 10~^ times smaller than the 
fundamental oscillation at Q.^. 

Stable solutions of eqs. (|2]) and (|3|) correspond to the 
upper and lower branches of a bistable hysteresis cycle. 
As the network analyzer only sweeps frequencies upward, 
curves in Fig. [2] (left) explore the upper branch until 
they suddenly fall down to the lower branch. To sweep 
frequency either up or down, we make use of a signal 
generator together with a spectrum analyzer in max-hold 
configuration (see setup h in Fig. [T]). Both sweeps are 
shown in Fig. [2] (right) and clearly exhibits an hysteresis 
cycle. 

The nonlinearity arises from additional stress induced 
by large displacements and we expect nonlocal effects as 
intrinsic elastic parameters such as Young's modulus of 
the membrane may be affected l3, [HI . We have thus 



investigated the influence of one mode (mode 1) on the 
natural frequency of another one (mode 2), using both 



the signal generator and the network analyzer as shown 
in Fig. [H 

An upward frequency sweep is first produced by the 
signal generator with a 20-dBm power to actuate the first 
mode in its nonlinear regime. In order to explore different 
points on the upper branch of the nonlinearity, the sweep 
is stopped at different frequencies l^p, represented by the 
arrows in Fig. [3] (left). While the drive of mode 1 remains 
at that frequency, we monitor the linear response of mode 
2 using the network analyzer with a 0-dBm actuation 
power. Central curves in Fig. [3l obtained for increasing 
actuation frequencies l^p, clearly show that the resonance 
frequency ^Ig of mode 2 is shifted as the displacement 
amplitude of mode 1 gets larger and larger. Considering 
that mode 2 has a quality factor of 6 000, the mode is 
actually displaced 170 times its width. We also explore 
the lower branch of the nonlinear response of mode 1 
using an initial downward frequency sweep (Fig. [3]right), 
thus retrieving the bistable behavior of mode 1 on the 
resonance- frequency shift of mode 2. 

As the effective elastic parameters of the membrane 
are altered by the nonlinearity, one may also expect a 
modification of its dynamics, which we have investigated 
using a ring-down technique. We have thus monitored 
the motion of the membrane in phase space, defining the 
slowly varying quadratures Xi and X2 as: 



Xp(t) = Xi(t) cos(ytQt) + X2{t) sm(ytQt). 



(4) 



Xi and X2 are readily accessible using the spectrum an- 
alyzer in IQ demodulation mode. Figure |4] shows the 
forced regime and the ring-down decay of the mechanical 
mode for different actuation levels: in the linear regime 
at its natural frequency VLq (left column), in the linear 
regime out of resonance (central column), and in the non- 
linear regime (right column). The membrane is actuated 
for times t < 0, with an upward frequency sweep to get 
the mode on the upper branch in the nonlinear case, then 
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FIG. 3. Intermodal coupling: actuation of mode 1 in its non- 
linear regime (left curve) induces a frequency shift of the res- 
onance of mode 2 (central curves) when the actuation fre- 
quency Qp of mode 1 is increased along the upper branch 
(arrows). Right curves display the bistable behavior of the 
resonance frequency of mode 2 when Qp is either swept up 
or down. The nonlinearity strength of mode 1 is comparable 
to the one experienced in curve e in Fig. O 
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FIG. 4. Time evolution in phase-space for the forced (blue) 
and free (red) regimes. Three cases are investigated, as 
sketched on top of each column, from left to right: linear 
regime actuated at resonance (5 dBm at 1057 kHz), out 
of resonance (9 dBm at 1056 kHz), and nonlinear regime 
(29 dBm at 1 058 kHz). Curves from top to bottom represent 
the phase-space trajectories and the amplitude and frequency 
evolutions. 



the actuation is stopped by a fast switch at time t = 0. 

The phase-space trajectories (top curves) first show the 
forced regime which appears as blue circles or as a single 
stationary point in the resonant case (left), and second 
the ring-down regime displayed as red spirals which re- 
duce to an almost straight line down to the origin when 
the motion is close to resonance. 

/^From these data we can infer the time evolutions of 
both the amplitude and frequency of motion (middle and 
bottom curves). In contrast to other systems such as 
graphene [l9|, the amplitude decay time is the same in 
the different actuation regimes: r = l/F 0.56 ms, cor- 
responding to a mechanical quality factor Q — 3 700. 
This tends to confirm that the damping in these res- 
onators does not depend on the effective Yo ung modulus 
and is indeed mostly due to surface effects |ll|. 

Yet the frequency behavior drastically changes be- 
tween the linear and nonlinear regimes. In the former 
case, even for a non-resonant actuation, the frequency in- 
stantaneously jumps back to its natural value when the 
actuation is released, as can also be inferred from the 
straight line trajectories in phase space. This result is 
actually expected for a linear harmonic oscillator, as its 
free evolution for t > only depends on the initial posi- 
tion and speed at t = 0, and not on the previous forced 



actuation frequency. In the nonlinear case, the frequency 
hops towards the top of the upper bistability branch and 
then slowly relaxes down to the natural frequency as the 
nonlinearity fades out. This can be understood from eq. 
([2]) which shows that for a given oscillation amplitude x^ , 
the preferred oscillation frequency is l^o (1 + ^/2) [1^. 
As the oscillation amplitude and then the nonlinear co- 
efficient s [eq. ([3|)] exponentially decrease with time, one 
thus expects the same behavior for the free oscillation 
frequency. 

To confirm this interpretation, we have studied the 
mechanical response of the mode to a weak probe ex- 
citation while it is simultaneously set in the nonlinear 
regime by a strong pump actuation. The setup is actu- 
ally similar to the one used for the two-mode actuation 
(Fig. [3]), except that both pump and probe are actuating 
the same mechanical mode, close to its natural frequency 
Qo/27r ~ 782 kHz for the membrane used here. Figure [5] 
(left) shows the spectra Xs[l^s] obtained when the probe 
frequency l^s is scanned around the pump frequency l^p, 
for different pump actuation powers. We clearly observe 
two symmetrical lorentzian peaks whose frequency differ- 
ence increases with the pump power. 

Using a spectrum analyzer instead of a network ana- 
lyzer, we actually checked that the two resonances simul- 
taneously appear: when actuated at a frequency l^s, the 
pump-probe combination leads to the concomitant gen- 
eration of a mechanical response at frequencies l^g and 
2Qp — l^s, symmetrically disposed on both sides of the 
pump frequency. 

This behavior can be understood from the Duffing 
equation ([1]), adding a second source term as{t) = 
2asCos(l]s^) related to the probe. The cubic term in 
Duffing equation is then responsible for the mixing of 
pump and probe signals, in a way similar to a x^^^ non- 
linearity in optics. Indeed, writing the total displacement 
in presence of the probe as the sum of the unperturbed 
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FIG. 5. Mechanical response of a mode pumped in the non- 
linear regime (frequency f^p) and probed by a weak actua- 
tion (frequency Qs). Spectra obtained with pump actuation 
powers of 10, 20, and 30 dBm (curves a to c) exhibit two res- 
onances symmetrically disposed around the pump frequency 
Qp , as demonstrated by curves d and e: dots are experimental 
values of the two resonance frequencies obtained for a pump 
actuation of 5 dBm (leading to a nonlinearity strength similar 
to curve c in Fig. [2]), solid curves are theoretical fits. 
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displacement Xp(t) due to the pump only [eq. (j2j)] and 
an additional displacement Xs(t), we obtain the following 
equation (assuming \xs\ <C |^p|): 

x,{t) + Tx,{t) 3/3x^(t)] x,{t) = a,{t). (5) 

As Xp{t) is oscillating at frequency l^p, the nonlinear term 
Xp(t) contains both a static term and a fast oscillation at 
2Qp. The former is responsible for a modification of the 
linear mechanical response at frequency whereas the 
latter is formally similar to phase conjugation in optics 
and generates a response at 2Qp — Qs- From ^ one 
gets the following coupled equations for the two Fourier 
components of x^ at frequencies and 2Qp — Q^'. 

C[n,] x,[n,]^enl x^,[2n^ - n,] = a,, (6) 
C[2% - n,] x^^[2% - ns]^enl x,[^,] = o, (7) 

where = + — — i^^o^ and we have as- 

sumed for simplicity Xp real. We thus obtain two phase- 
conjugate lorentzian responses with resonance frequen- 
cies: 

(7s+ = no{l^e), Qs- = 21^p - a+. (8) 

Note the difference between and the resonance 

frequency Qo{l-\-e/2) in response to a single actua- 
tion [eq. (j2j)], actually equivalent to the difference ob- 
tained between self- and cross-phase modulations in Kerr 
media [20|. 

Dots in Fig. [5] (right) show frequencies and Qs- 
deduced from spectra similar to the ones displayed on 
Fig. [5] (left), as a function of the pump frequency Qp. 
Curves d and e are theoretical predictions obtained using 
eqs. dl]) with s estimated from the experimental nonlin- 
ear behavior [eq. (j2j) and Fig. [2]. At low pump frequency 
r^p ^ r^O: the membrane is almost in a linear regime 
and we only observe one peak close to the natural fre- 
quency flo' At higher frequency, the two resonances are 
symmetrically located around ftp (dashed hue), in good 
agreement with theoretical predictions. 

We have demonstrated nonlinear effects in the dynam- 
ics of PhC nanomembranes: bist ability, intermodal tun- 
ing and the generation of phase conjugate modes. Pump- 
probe experiments have in particular allowed us to un- 
derline the non-local character of the Young modulus 
modification. Such nonlinear behaviours are intrinsic to 
nano- and microscale systems Q, and have dramatic con- 
sequences not only in the fields of optomechanics and 
nonlinear dynamics, but in quantum optics as well, as in- 
termodal coupling can be used to perform QND measure- 
ments of a signal mode by monitoring a meter one 
[23| . Optomechanical systems can also be used for classi- 
cal or quantum information processing, either by taking 
advantage of their long mechanical coherence times to 
store information in a mechanical excitation [2^, or by 
using their nonlinear character for non volatile memories 
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